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Using Mandelstam's approximation to the gluon Dyson-Schwinger equation we 
calculate the gluon self-energy in a renormalisation group invariant fashion. We 
obtain a non-perturbative /3 function. The scaling behaviour near the ultraviolet 
stable fixed point is in good agreement with perturbative QCD. No further fixed 
point for positive values of the coupling is found: as increases without bound in 
the infrared. 



The infrared behaviour of the running coupling in strong interactions, as, 
is of great interest, since it may provide an understanding of confinement. Its 
study is an intrinsically non-perturbative problem. One suitable framework to 
address this problem is provided by the Dyson-Schwinger equations of QCD. 
Studies of this infinite tower of equations rely on specific truncation schemes. 
Here we will focus on an approximation to the gluon Dyson-Schwinger equation 
originally proposed by MandelstamH This yields a simplified equation for the 
inverse gluon propagator in Euclidean momentum space, 



D-^{k)=D^{k) + Jg^Jd*qrr(k,q-k,-q) (1) 

xD^(q)D^(k - q)T^(q, k - q, -k) , 

where Do and Tq are the bare gluon propagator and the bare three-gluon 
vertex. The use of a second bare vertex in the gluon loop is combined with 
one of the gluon propagators being bare. This is to account for some of the 
dressing of the vertex as entailed by its Slavnov- Taylor identityu 
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In a manifestly gauge invariant formulation the Dyson-Schwinger equation 
for the inverse gluon propagator in the covariant gauge would be transverse 
without further adjustments. This may be violated due to the neglect of ghosts, 
the violation of Slavnov-Taylor identities and also due to a regularisation that 
does not preserve the residual local invariance under transformations generated 
by harmonic gauge functions (d 2 A(x) = 0). The latter is the case for an 
0(4) invariant Euclidean cutoff A, which we will use to regularise eq. (0). 
Contracting eq. (0) with the transversal projector, 



P ^ k ) = - Z— . (2) 



k^k h 
~~k 2 

leads to an equation for the gluon renormalisatiori function G(k 2 ) in Landau 
gauge, which after the angular integrations read£la 

G^)=l +I^lHo fc V(f£-¥f^§)G( 9 2 ) (3) 

I gn 1 f A2 J„2 / 7fc* _ 25 l£ _ I) r(n 2 \ 
^IQir 2 k 2 Jk 2 ™ \8 9 4 4 q 2 2 J U W I ■ 

However, the above equation contains a quadratically ultraviolet divergent 
term, which has to be subtracted by a suitable counter term. Generally, 
quadratic ultraviolet divergences can occur only in the part of the inverse gluon 
propagator proportional to 8^ . Therefore, that part cannot be unambiguously 
determined, it depends on the routing of the momenta. This is due to the var- 
ious violations of gauge invariance mentioned above. The unambiguous term 
proportional to k^k v can be obtained by contracting (|l|) with 

(k) =5^-4-^- . (4) 
k 

In this case, upon angle integration instead of (^) one obtains! 

' 1 +li^/o fc2 ^(l^-¥^-l)^ 2 ) (5) 



G(k 2 



The logarithmic ultraviolet divergences in eqs. 

d) and (§) can be removed 
by multiplicative renormalisation. Introducing renormalised gluon propagator 
and coupling the renormalisation constants Z% and Z g are defined by D — > Z3D 
and go — Z g g. The renormalised Dyson-Schwinger equation for the gluon 
propagator in Mandelstam approximation then reads, 

n -1 



G(k 2 ) = 



(G) 
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Figure 1: The gluon renormalisation function gG(k 2 ) = 8ir<j/k 2 + gF(k 2 ) for eq. Q The 
dashed lines show gF(k 2 ) for eqs. Nand ra. 



with obvious definitions of Iq as a functional of G for the respective cases (c/., 
eqs. |^ and ||). In Mandelstam approximation one has Z g Z% = la We adopt 
a momentum subtraction scheme requiring the gluon self-energy to vanish at 
the renormalisation point /i: G(fi 2 ) = 1. 

The behaviour of G(k 2 ) for k 2 — > can be summarized as follows: 



gG(k 2 ) = "— +9 F{k 2 ) , ^(fc 2 )=a 00 (fc 2 /A 2 ) T + ... . (7) 
Here, a is the string tension, and the subleading term is determined bj 



/31 

y=\ Irs 1.273, and a 00 m 0.29421 

V 6 

for Mandelstam's original equation (^). For eq. (g), we obtain, 

2 / 1 ( 1099 \\ 13 

7 = - V229cos - arccos == 1.271 , and a 00 ~ 0.29446 . 
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The numerical solutions to eqs. (|^) and (|^) for the renormalisation group in- 
variant products gG and gF as functions of fc 2 /A 2 vIOM are shown in fig. [j]. As 
there are no qualitative and only little quantitative differences between those 
two solutions, we concentrate on the discussion of the solution to eq. (^|) in the 
following. 
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The non-perturbative and renormalisation group invariant result relates 
the string tension a to the QCD scale (the only parameter) by a = 2A 2 . Fixing 
a to its phenomenological value (it may be extracted from heavy quarkonium 
spectra or the slope of Regge trajectories) gives a value for A of about 600 MeV. 
Alternatively, from as{m 2 z ) ~ 0.108 we find A ~ 640MeV. For a momentum 
subtraction scheme with Nf = this is the correct order of magnitude. 

The scaling behaviour of the solution near the ultraviolet fixed point is 
determined by the coefficients f3o — 14, j3\ = 70/3 and 7J4 = 7 which are 
reasonably close to the perturbative values for Nf = 0, i.e., /3q = 11, j3\ = 51 
and 7^ = 13/2. 

The running coupling, which is obtained for arbitrary scales from the renor- 
malisation condition, resembles the two-loop perturbative form asymptotically. 
It is shown together with the corresponding Callan— Symanzik (3 function in 
fig. |[ The running coupling increases without bound in the infrared, no fur- 
ther fixed point exists. 
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